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I. INTRODUCTION 



Massive color-octet vector bosons are predicted in a variety of models, including axigluon models [HE], topcolor 
models [3 -6 , technicolor models with colored technifermions [7 , flavor-universal [8 , 9 and chiral [10] coloron models, 
and extra-dimensional models with KK gluons [Tl]|12]. These states have also recently been considered as a potential 
source [I3l of the top-quark forward-backward asymmetry observed by the CDF collaboration [151 [IB| {j Recent 
searches for resonances in the dijet mass spectrum at the LHC imply that the lower bound on such a boson is now 
2-3 TeV [2OH2310 If there are color-octet vector bosons associated with the electroweak symmetry breaking sector, 
as suggested by several of the models discussed above, their presence should be uncovered by the LHC in the future. 

In this paper, we report the first complete calculatiorj^ of QCD corrections to the production of a massive color- 
octet vector boson. We will refer to these massive color-octet vector states generically as "colorons." We treat the 
coloron as an asymptotic state in our calculations, employing the narrow width approximation. Our next-to-leading- 
order (NLO) calculation includes both virtual corrections as well as corrections arising from the emission of gluons 
and light quarks, and we demonstrate the reduction in factorization-scale dependence relative to the leading-order 
(LO) approximation used in previous hadron collider studies. 

The QCD NLO calculation of coloron production reported here differs substantially from the classic computation 
of the QCD NLO corrections to Drell-Yan production [26], because the final state is colored. In particular, Drell-Yan 
production involves the coupling of the light quarks to a conserved (or, in the case of W- or Z-mediated processes, 
conserved up to quark masses) current. Hence, in computing the NLO corrections to Drell-Yan processes, the 
current conservation Ward identity insures a cancellation between the UV divergences arising from virtual quark wave 
function and vertex corrections. These cancellations do not occur in the calculation of the NLO corrections to coloron 
productio n, b ecause of vertex corrections involving the 3-point non-Abelian colored-boson vertices. As we describe 
in Section IV, we use the "pinch technique" [27 to divide the problematic non-Abelian vertex corrections into two 
pieces - a "pinched" piece whose UV divergence contributes to the renormalization of the coloron wavefunction (and, 
ultimately, a renormalization of the coloron coupling) and an "unpinched" part whose UV divergence (when combined 
with an Abelian vertex correction) cancels against the UV divergences in quark wavefunction renormalization. As 
we show, once the UV divergences are properly accounted for, the IR divergences cancel in the usual way: the IR 
divergences arising from real quark or gluon emission cancel against the IR divergences in the virtual corrections, and 
the IR divergences arising from collinear quarks or gluons in the initial state are absorbed in the properly defined 
parton distribution functions (PDFs). 

We compute the gauge-, quark-, and self-couplings of the coloron from a theory with an extended SU{3)ic x 
SU{3)2c SU{3)c gauge structure, where SU{3)c is identified with QCD. The calculation yields the minimal 
coupling of gluons to colorons, and allows for the most general couplings of quarks to colorons. The cancellation of 
UV divergences described above, however, occurs only when the 3-coloron coupling has the strength that arises from 
the dimension-four gauge-kinetic energy terms of the extended SU{3)ic x SU{3)2c gauge structure. Our computation 
applies directly to any theory with this structure, i.e. to massive color-octet vector bosons in axigluon, topcolor, and 
coloron models. In general, the triple coupling of KK gluons in extra-dimensional models, or of colored technivector 
mesons in technicolor models, will not follow this pattern. However our results apply approximately to these cases as 
well, to the extent that the SU{3)ic x SU{3)2c model is a good low-energy effective theory for the extra dimensional 
model (a "two-site" approximation in the language of deconstruction [28, 29 ) or for the technicolor theory (a hidden 
local symmetry approximation for the effective technivector meson sector [30, 31 )|^ 

This paper is structured as follows. In Sec. [llj we introduce the formalism of a minimal vector coloron theory. 



deriving all the Feynman rules, and setting the stage for the subsequent calculations. In Sec. HI we review the 



leading order computations of the amplitude and cross section for coloron production due to qq pair annihilation. 
Sec. |IV| describes in detail the one loop virtual corrections to the qq pair annihilation process, elaborating on the 
contributions from the quark self-energy, coloron-coloron and gluon-coloron mixed vacuum polarization amplitudes, 
and the vertex corrections. We employ the pinch technique [27 , described above, in order to consistently treat 
the UV divergences, and obtain a gauge- invariant, mutually-independent set of counterterms. The one loop cross 
section is constructed, and the IR singularities of the virtual correction properly extracted. In Sec. [Vj we consider 



the real emission processes, consisting of real (soft and collinear) gluon and (collinear) quark emission. In Sec. VI 



we put all the pieces together, exhibiting the explicit cancellation of the IR divergences among the real and virtual 



^ Note, however, that the observation of a top-quark forward-backward asymmetry is not confirmed by results of the DO collaboration 
[TTIIII]- Furthermore, if the observed top-quark forward-backward asymmetry is confirmed, explaining this using color-octet vector 
bosons is problematic given the tight constraints on fiavor-changing neutral-currents |19| . 

^ At least for the fermion charge assignments considered, and in the case where the resonance is narrow compared to the djiet mass 
resolution of the detector 

^ As this work was being completed, a computation of the NLO virtual corrections of top-quark pair production via a heavy color-octet 
vector boson has been reported in 24 . That work is complementary to ours in that it does not employ the narrow width approximation 
for the color-octet boson, but neither does it include real gluon or quark emission. After this work was submitted for publication, real 
emission has also been considered by those authors 25^. 

Arbitary three- and four-point coloron self-couplings can be incorporated in the SU{3)ic x SU{3)2c by adding 0{p'^) terms in the of 
effective chiral Lagrangian of Eq. and deviations in these couplings are therefore of 0{M^ / where A is the cutoff of the effective 
coloron theory. The 3- and 4-point self-couplings, however, are neither relevant to the leading-order qq nor to the IR divergent NLO 
coloron production contributions, and therefore numericallv insignificant. 
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corrections, and demonstrate the renormalization of the quark and gluon PDFs. We give a finite expression for the 
NLO-corrected production cross section. Finahy, in Sec. |VII| we plot the cross section, demonstrate that the QCD 
NLO corrections are as large as 30%, and show that the residual factorization-scale dependence is at the 2% level. 
We also calculate the i^-factor and the pt spectrum for coloron production, since these are valuable for comparison 
with experiment. 

An appendix contains all the Feynman rules of the theory. 



II. A MINIMAL THEORY FOR SPIN-ONE COLORONS 

In this section, we introduce coloron^ as the massive color-octet bosons arising when an extended 6'/7(3)ic x 
SU (3) 2c gauge symmetry is spontaneously broken by a non-linear sigma model field to its diagonal subgroup, SU (3)c, 
which we identify with QCD. The symmetry breaking results in a low-energy spectrum that includes both a massless 
spin-one color octet of gauge bosons, the gluons, and a massive spin-one color octet of gauge bosons, the colorons. 

In detail, we replace the QCD Lagrangian with 

>^color — ~ ~^G\^uG^ — ~G2^v'G2 H ~ TrD^SD^S^ + >^gauge— fixing ~^ >^ghost ~^ ^quark • (1) 

Here S is the nonlinear sigma field breaking SU{3)ic x SU{3)2c to SU{3)c, 

S = expf^^j , a = l,...,8, (2) 

where tt^ are the Nambu-Goldstone bosons "eaten" by the coloron, / is the corresponding "decay-constant", and 
are the Gell-Mann matrices, normalized as Tr t^t^ = The E field transforms as the bi- fundamental of 

SU{3)i, X SU{3)2c, 

H , Ui = exp {ia'^t^) , (3) 

where the are the parameters of the SU{3)ic transformations. This leads to the covariant derivative 

I)^S = 5^S-z^,,Gy^S + i^,,SGy^ , (4) 

where ^g. is the gauge couphng of the SU{3)ic gauge group. Up to a total divergence, the quadratic terms in the 
Lagrangian are 

+ £(2) ,£(2) ,£(2) /5^ 

' gauge — fixing ' ghost ' quark ' \ / 

where a sum over z = 1, 2 in the gauge kinetic terms is implied. 

The gauge-Goldstone mixing term can be removed, up to a total divergence, by choosing the gauge-fixing La- 
grangian to be 

>^gauge-fixing = ) ' 

where the gauge-fixing functions are 



The Faddeev-Popov ghost Lagrangian is obtained by taking the functional determinant of ST^/6a^. This leads to 
/-» — ( 



.91 f 



- gs, r'^jGl) - - Si2) {Si, - 52,) 5^' + 0(7r)J c] , (8) 



where /"^^ are the SU{3) structure constants, and a sum over i,j = 1,2 is implied. Notice that we have included 
only the inhomogeneous terms in the transformation of the "eaten" Goldstone boson, whence the unspecified 0{7r) 



^ Colorons can in principle be introduced as matter fields in the adjoint of SU(3)c- This approach, however, would lead to an early 
violation of tree-level unitarity, as the scattering amplitude of longitudinally polarized massive spin-one bosons can grow, by power 
counting, like E"^, where E is the center-of-mass (CM) energy. The only way to avoid this is to "promote" the coloron to the status of 
gauge field of a spontaneously broken gauge theory: then the special relation between trilinear and quartic gauge couplings will lead 
to an exact cancellation of the terms growing like E'^, as happens in the standard electroweak theory. 
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terms in the ghost Lagrangian, which are unnecessary for our computation. Up to a total divergence, the quadratic 
Lagrangian now reads 



C 



(2) ^ 
color 



Gl 



Sijg'^'d^ - Sij 1 



g^Q. + ^ {da - Si2) (<5y - S2j) 



9lP 



,(2) 



r°- -I- n 
3 '-'quark • 



(9) 



Aside from a factor of the gauge-fixing parameter ^, the gauge and ghost fields share the same mass matrix, as 
expected. This is diagonalized by 



1^ \ - R 



1 - 



R 



where 



R 



cos e/c — sm ( 
sin On cos 0, 



smf 



9s^ 



(10) 



(11) 



In Eq. ( 10 ) is the gluon field and is the coloron field, whereas Cg. and are the corresponding ghost fields. 
Inserting these expressions in Eq. ((9| gives, for the coloron mass. 



Mc 



\/9i, + 9% f _ 9s f 



(12) 



where Qs is the S'/7(3)c coupling. 



91 9l, 



7^ 

^S2 



(13) 



The gluon ghost is massless, whereas both the coloron ghost and the eaten Goldstone boson have mass ^/^Mc^ The 
interaction vertices and the corresponding Feynman rules can be found in Appendix [Aj 

We will leave the quark charge assignments under SU{3)ic x SU{3)2c arbitrary, for greater generality. In the mass 
eigenstate basis we write 



/:quark = Q^i [0 - i9s0"t^ - 10" {9lPl + 9rPr)\ q^ , 
where Pl and Pr are the helicity projection operators. 



1-75 



1+75 



(14) 



(15) 



2 ' 2 ' 

and Hs a flavor inde:?^ The coupling to the gluon is dictated by charge universality, whereas the ql and qr couplings 
to the coloron depend on the original charge assignments of the quarks. For example, if both left-handed and right- 
handed quarks are only charged under S'/7(3)ic, then 9l = 9r = — ^s^an^c^ while the axigluon [1, 2 corresponds 
to = —9r = 9s (^-e- = 7r/4). In general, gL and gR can each take on the values —^5 tan 6>c or gs^otOc in any 
specific modelQ 



9l,9r ^ {-^s tan l9c,^s cot6>c} 



(16) 



III. LO COLORON PRODUCTION 

The dominant channel for coloron production at a hadron collider is given by the tree-level diagram of Fig. [l] 
in which a qq pair annihilates into a coloron. The tree-level diagram with gluon-gluon fusion into a coloron does 
not exist in the Lagrangian of Eq. ([T]): in general there are no dimension-four terms with two gauge bosons of an 
unbroken symmetry and a spin-one field charged under the same symmetry. We use the narrow width approximation 



^ Here we work in the broken electroweak phase, and only employ fermion mass eigenstates. 

^ It is possible to generalize this setup to non-universal charge assignments: in this case flavor-diagonal chiral couplings to the coloron 
would depend on a generation index. Flavor-changing couplings are strongly constrained [19| . 
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a, V 




FIG. 1. Tree- level contribution to coloron production. The coloron field, C^, is represented by the zigzag line. 



for the coloron, take the quarks to be on-shell, and set their masses to zero: this is certainly a good approximation, 
as current experimental bounds [20^, '22', ^23^ constrain the coloron mass to be in the TeV range. 
The leading order (LO) amplitude corresponding to the diagram of Fig. [l] is 

i>(iy_^c = 9s v^iP) {tlPl + trPr) e u\p) ef*{r) , (17) 
where the superscripts r and s denote quark spin projections, A is the coloron polarization, and 

tl = ^ , rR = ^ , vl^tr e {- tan 6>c, cot 6>c} . (18) 
9s 9s 

In d = 2(2 — e) dimensions the squared amplitude averaged over initial spins and colors, and summed over final 
polarization states, is 

^ ^ ^ ^ ^ ^ spin & color ^ ^ 

where dim(r) = 3 and C2(r) =4/3 are respectively the dimension and Casimir of the fundamental representation of 
5'[/(3), and s = {p -\- p)'^ = 2p - p is the partonic center-of-mass CM energy. This gives the LO cross section [2 for 
qq C, 



-Stc = i \M^^c\' Sil - X) = ^^^±^ Sil - X) , (20) 



where as = ^g/47r. 



and 



X.^. (22) 



The full LO cross section for pp ^ C is given by the convolution of the LO partonic cross section d-^^^c with the 
parton distribution functions (PDFs) for the quarks within the protons, and a sum over all quark flavors. 



J dXi J dX2Y^ [fq{Xl)fq{x2)^fq{xi)fq{x2) Cff^^^ , (23) 



where fq{x) is the PDF of parton and x the momentum fraction of the corresponding parton. Taking the collision 
axis to be the 3- axis, the four- moment a of the partons are: 

p= ^(xi,0,0,xi) , p= ^(X2,0,0,-X2) , (24) 
where s is the CM energy of the colliding hadrons. This gives 

s = XiX2 s , x= ^ • (25) 

S Xi X2 
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FIG. 2. Structure qq ^ C amplitude, to all orders in perturbation theory. Direct coloron production is illustrated on the 
left, while production via mixing with the gluon is shown on the right. The gluon field is, as usual, represented by the coiling 
line; the coloron field is represented by the zigzag line. 




P P p p 

FIG. 3. Quark self-energy diagrams at one loop. Particle notation as defined in Fig. [2] 



IV. NLO COLORON PRODUCTION: VIRTUAL CORRECTIONS 



In this section we compute the next-to-leading order (NLO) virtual QCD corrections to the qq ^ C amplitude. 
These include one- loop wave- function and vertex corrections, which we choose to compute in 't Hooft-Feynman 
gauge, ^ = 1. The non-Abelian vertex corrections are computed by employing the pinch technique: this allows us to 
obtain QED-like Ward identities, and absorb all UV infinities in the renormalization of the gauge field propagators. 
After inclusion of the counterterms, the virtual corrections are UV- finite, yet IR infinite. In Sec. [V| we show that 
the IR divergences cancel once the real corrections, corresponding to the emission of soft and collinear gluons and 
quarks, are included in the calculation of the inclusive production cross section. Our loop integrals are computed in 
dimensional regularization, with d = 2{2 — e) dimensions. We first regulate the IR divergences by giving the gluon 
a small mass {rrig ^ in this way all infinities are in the UV, and regularization requires e > 0. After all of the 
UV infinities are removed, by cancellation and inclusion of the counterterms, we let the gluon mass approach zero. 
This will make the virtual corrections IR divergent, with the infinities being regulated by taking e < 0. 

Since the quark couplings to the coloron are chiral, in general, we need a prescription for treating 75 in d 7^ 4. 
Here we take 75 to always anticommute with 7^. Choosing an alternative prescription, such as 't Hooft-Veltman in 
which 75 anticommutes with 7^ for /i = 0, 1, 2, 3 and commutes for other values of /i, would lead to a cross section 
for qq ^ C which differs from ours by only a finite renormalization of the coupling(s). 

The general structure of the qq ^ C amplitude, illustrated in Fig. |2j is 

Z,u%p) ef*{r) , (26) 

where T^qc i^qqc) one-particle- irreducible (IPI) quark- quark- coloron (quark-quark-gluon) vertex and Hoc is 

the coefficient of g^^ in the gluon-coloron vacuum polarization mixing amplitude (VPA). The factors Zq and Zc are, 
respectively, the residues of the full quark and coloron propagators at the mass pole; they are obtained from the 
quark self-energy amplitude, S(^), and the coefficient of g^^ in the coloron-coloron VPA, IIcc{q'^)^ as follows: 



iM 



qq^C 



v^{p) i 



^ qqC ^ qqG 



FIG. 4. One- loop Abelian vertex correction to the qq ^ C amplitude. Particle notation is as defined in Fig. [2] 



where the prime denotes a derivative with respect to the argument, and M^phys is the coloron's physical mass. To 



lowest order, Zq = 1, Zc = 1, ^gc = 0, and iT^qc ~ 7^ {'^lPl + trPr) t^; inserting these in Eq. (26) recovers the 



tree- level amplitude of Eq. (17). 



A. Quark Self-Energy 



The NLO quark self-energy correction to the qq ^ C amplitude is found, from Eqs. (26) and (27), to be 



where 



iQ = v-{p) iY {gLPL + gnPR) SZ, u'ip) *(r) , 



5Zg = S'(0) . 



(28) 
(29) 



At one loop, the S(|i) amplitude is given by the diagrams of Fig. |3] These lead to the expression 



S(i^) = -p 



5,2C2(r)2(l-e)r(e) f\,^ ^,\(^^Y . 



167r2 



dx (1 — x) 



V A, 



Gq 



^Cq 



(rlPL+rlPn) 



(30) 
(31) 



where r(e) is the Euler Gamma-function evaluated at infinitesimal e, and 

= (1 — x)m^g — x(l — x)p^ — if] , Acq = (1 — x)Mq — x{l — x)p^ — irj . 

The parameter ja is the mass scale introduced by the loop integral in d dimensions, and r] is the positive infinites- 
imal parameter giving the appropriate prescription for computing the integral in momentum space. As previously 
anticipated, we have introduced a small gluon mass, m^, in order to regulate the IR divergences and isolate the UV 
infinities: with rrig ^ 0, ^{p) and S'(^) contain only UV divergences. Inserting Eq. (30) in Eq. (29) gives 



SZn 



167r2 



dx (1 — x) 







47r/i^ 



The amplitude of Eq. ( |28[ ) becomes 



^ 2C2(r)(l-e)r(e) / dx 
47r Jo Jo 



dy 



(1 - x)ml 
47r/i^ 



47r/i^ 



(1 - x)M^ 



(1 — x)m^ 



IT] 



iM 



(0) 



{tIPl^tIPr) 



.(32) 



qq^C^ \ {l-x)Ml-ir] 



where M^^^^^j is given 



by Eq. (|l7|). 



and 



(33) 



(34) 



For later convenience we have traded the 1 — x factor, in Eq. ( |32[ ), for an integral over dy: this will allow us to directly 
add the self-energy correction to the vertex correction and explicitly show the cancellation of the UV divergences. 



B. Abelian Vertex Corrections 



The one-loop Abelian vertex correction to the qq ^ C amplitude is given by the diagrams of Fig. |4l These lead 
to the amplitude 



a. 



^eiian = 7z[2C'2(r)-C2(G)]r(l + e)^ dxj^ 



47r 

(1-6)2 



{xye - (1 - x){l - y)) 



^Cqq 



/ 4V 

K^Cqq 



iM 



{xye - (1 -x){l -y)) 



/(O) 
qq^C 



^Gqq 



/ 47r/i 

K^Gqq 



2 \ e 



iM 



(0) 

qq^C 



(35) 



FIG. 5. One-loop non-Abelian vertex correction to the qq ^ C amplitude. Particle notation is as defined in Fig. [2] 
three-gauge-boson vertex in these diagrams is a full non-Abelian vertex F^^^ in Eq. (37). 



Each 



where C2{G) = 3 is the Casimir of the adjoint representation, and 

^Gqq = (l-x- y)ml - xys - ir] , 

^Cqq = {l-x- y)M^ - xys - if] . (36) 
Once again, we have included a small gluon mass rUg in order to regulate the IR divergences. 

C. Non-Abelian Vertex Corrections a la Pinch- Technique: Unpinched Diagrams 

The non-Abelian vertex corrections are given by the diagrams of Fig. [5] When added to the overall Abelian vertex 



correction, Eq. (35), these give the one-loop total vertex correction to qq C . Unlike in QED, the UV divergences in 
the vertex correction do not cancel the UV divergences arising from the self-energy amplitudes. The reason for this is 
that the QED Ward identity d^j^ = is now replaced by its non-Abelian counterpart D^j^ = 0, which does not imply 
the equality of vertex and quark-wavefunction renormalization constants. It is possible, though, to recover QED-like 
Ward identities for the currents by employing the pinch technique. This consists of breaking up the gauge-boson 
internal momenta of a Feynman diagram into "pinching" and "non-pinching" pieces. The pinching momenta are 
those which cancel some internal propagators, leading to a simpler diagram with the external-momentum structure 
of a propagator. The non-pinching momenta will instead give overall amplitudes satisfying QED-like Ward identities. 
A formal proof of these statements, for an arbitrary non-Abelian gauge theory, can be found in the review of Ref. [ 27] 
(and references therein). 

In our vertex computation the pinch technique works as follows. The non-Abelian vertex structure in each of the 
diagrams in Fig. [5] is 

F^"^(/c,p, p) = g^^{-2p -p^ky^ g^P{p -p- 2kY + gP^{k + p + 2pY . (37) 
We can break this into two parts, 

T^^^''{k,p,p) = T%^''{k,p,p) + T^;''{k,p,p) , (38) 

where 

r''/%k,p,p) = -2g>^''ip + pr + 2gP>^{p + pr + g''P{p-p-2kr , (39) 
T'ir''{k,p,p)=g>'''{p + ky+gP>'{k-pr ■ (40) 

Unlike r''''P{k,p,p), the T^p''{k,p,p) vertex satisfies a QED-like Ward identity for tfie gC ^ C and CC ^ C 
amplitudes, 

{p + P)^T%'''{k,p,p)=g'^''[{p-kf-{p + kf] . (41) 

As shown below, when {k^p^p) is used to compute the integral in momentum space (instead of T^^^{k^p^p))^ its 
UV divergences, added to the UV divergences of the Abelian vertex corrections, exactly cancel the UV divergences 
of the quark self-energy amplitudes. As mentioned above, this occurs because a QED-like Ward identity for qq ^ C 
holds, as one can prove by using the QED-like Ward identity for the gC C and CC ^ C amplitudes given in 



Eq. (41 ). The three diagrams which correspond to using {k^p^p) instead of r^^^(/c,p,p) are symbolically denoted 
with ablack disk over the non-Abelian vertex, and are shown in Fig. |6] These lead to the following contribution to 
the qq ^ C amplitude: 






FIG. 6. Non-Abelian un-pinched vertex-correction diagrams for the qq ^ C amplitude at one loop. Particle notation is as 
defined in Fig. [2] The black disk indicates that each three-point gauge boson vertex in these diagrams has been replaced by 
the non-pinched portion F^^^, as described in Eqs. (38) and (39). 



iVr 



non— Abelian 



1 



C2(G)F(l + e) / dx 





e 

1 - e 



(x^y) 




U 



GCq 



^CCqJ \^CCqJ 



iM 



(0) 

qq^C 



^CCqJ \^CCq 



(42) 



where 



^GCq = xml + yMc - xys - ir] 
^CGq = xM^ + yml - xys - if] 
^ccq = + y)Ml. - xys - if] . 



(43) 



In order to obtain Eq. (42) we have used the equations of motion for the external spinors, together with the relations 

2cot2l9rL 



-1 



2 cot 2(9 tr = -1 



' R ^ 



(44) 



which are true for any charge assignment of the quarks. As anticipated, iQ + ^VAbeiian + ^Kon-Abeiian is free of UV 
divergences, as manifestly shown by adding together Eqs. (33), (35) and (42). This part of the amplitude is however 
IR divergent in the limit of zero gluon mass. Setting rUg = and e < gives 



iQ + ^VAbelian + ^Kon — Abelian 



47r 



C2{r) 



3 + 2i 



C2{G)- 



finite . 



(45) 



Of course we still need to include the contribution from T^^{k^p^p) (of Eq. (40)) in the full non- Abelian vertex 
correction. This contains the pinching momenta: the action of p and p on the external spinors gives zero, and 
the remaining piece cancels the internal fermion propagator in the diagram. Thus the internal fermion line in each 
diagram is "pinched" away, leaving an effective diagram with a 4-point coupling between fermions and gauge bosons 
as shown in Fig. [7| The UV divergences of the pinched diagrams have the same group- and momentum-structure as 
those of the VPAs, and can be absorbed in the counterterms for the gauge field propagators. In order to see this 
clearly, we will now consider the form of the "true" propagator corrections to the qq ^ C amplitude in the following 
subsection. 



D. Form of the Vacuum Polarization Amplitudes 



The NLO corrections to the qq ^ C amplitude due to the VPAs are found, from Eqs. (26) and (27), to have the 
form 



where 



- ^J^lqq^c 2 ^-^^^qq^C - 



(46) 



SZc = U^cc{Mh) , 



(47) 
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and 



iMl'Ec = 5« h^'t'^ ^%p) ■ 



(48) 



In order to obtain the second term of Eq. (46), we have replaced F^^^ with its LO component ij^t^. Notice also 
that at this order we can swap M^^^iys 

At one loop IIcc{q'^) is given by the diagrams of Figs, [s] and 10, in which the gluon ghost is represented by 
dotted lines, the coloron ghost by a sequence of filled circles, and the eaten Goldstone bosons are represented by 
dashed lines. There are poles at d = 2 proportional to both q'^ and M^. The latter correspond to quadratic 
divergences (renormalizing the coloron mass scale /), whereas the former can only be logarithmic by dimensional 
analysis (renormalizing the coloron field) |^ The momentum dependent part of the full coloron-coloron VPA is not 
transverse, as the coefficients of the q'^ and q^q^ terms are different. However we have explicitly verified that the 
infinite part is transverse: this is necessary, because the corresponding Lagrangian counterterms are transverse. For 
small values of e we obtain 



Utt 



2 



ncc(g')^^" + q^q"" terms = C2(G) 
^ (^l-x{4-3x)^E -x{l- 




■2xf 



^GC 



2x E^3-5x 







+4cot^(2i9c)C2(G) 

x{l — X 



dx< 



^cc 



Jo 



4 
dx 



E 



x{l — X 



^ (l^4x{l-x)^E^{l-2xy 
(1 - 2x)2 



^cc 



E q^'q'' 



^CC 



g'"'Ml. 



2x{l - x)\e {g^''' q^ -q^'q'') 



(49) 



where our results depend only on the coefficient of the quantity Nf is the number of quark ffavors in the loop 
(see Fig. [10|, 



7 + log 47r 



e 

and 7 is the Euler-Mascheroni constant. The A functions in Eq. (|49| are 

Agc = xml^(l- x)Ml - x(l - x)q^ 
Acc = M^- x{l - x)q^ - if] , 
—x{l — x)q^ 



(50) 



IT] , 



^CC 

A 



qq 



IT] . 



(51) 



Notice that the coloron-coloron VPA of Eq. ( |49[ ) is not IR divergent in the limit rrig 0, since there are no 
contributions with only massless (gluon) states. However what enters in Eq. (46) is the derivative of Hcc (see Eq. 
(47)), which is IR divergent in the limit rrig 0. 



FIG. 7. Non-Abelian pinched vertex-correction diagrams for the qq ^ C amplitude at one loop. Particle notation is as defined 
in Fig. |2] 



^ This situation parallels the renormalization of the electroweak chiral lagrangian [321 133) . 
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FIG. 8. Coloron-coloron vacuum polarization amplitude at one loop. A gluon field is, as usual, represented by a coiling line; 
a coloron field is represented by a zigzag line. The coloron ghost is represented by a sequence of filled circles, and the eaten 
Goldstone bosons are represented by dashed lines. 
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FIG. 9. Gluon-coloron mixing amplitude at one loop. Particle notation is as defined in Fig. [S] 



The momentum dependent part of the gluon-coloron mixing amplitude (Figs. [9] and 10) is found to be transverse, 
both in the infinite and the finite parts. For small values of e we find 



(^) 'nGc(^')^^^ + ^^^^ terms 

= 2cot(2i9c)C2(G) / dx 
Jo 

+ {VL + rR)Nf f 
Jo 



(52) 



.2 \ ^ 



^CC 



dx 



There are no potential IR divergences hidden in Hoc- 



(^^^ 5x{l - x)^E ^ {1 - 2xf 



{q^g^^ - q^q^) ^ E M\ 



E. Non-Abelian Vertex Corrections a la Pinch- Technique: Pinched Diagrams 

The pinched diagrams of Fig. [7| are obtained from the diagrams of Fig. |5]by replacing the full non-Abelian vertex 
momentum structure T^^^{k^p^p) from Eq. (37), with {k^p^p) from Eq. ([4Q|). This leads to the amplitude 



pinched 



47r 
47r 



C2(G) / dx 





2('^y + 4cot2(2^,) ^ ^ 



2 \ e 



^CC 



EM 



(0) 

qq^C 



2 cot (2(9, 



Jo 



dx 



2 \ e 



^CC 



E M 



//(O) 

qq^C ' 



(53) 



the form of a VPA correction, like that in Eq. (46). In fact we can write 



where we have used Eq. (44) to rewrite the fermion couplings in terms of Oc- This contribution to the amplitude has 

(54) 



^^pinched = 



(0) KAMcl^-j^m ^GC{S 



qq^C- 



qq^C 



where 



Utt 



) 'ncc(g') = C2(G)^'dx(^^) A{q^-Ml)E 



+ 4cot^(2(9c)C2(G) / dx 



/' 

Jo 



^cc 



2{q^ - Ml)E 



and 



(55) 



(56) 
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FIG. 10. Fermion contributions to coloron-coloron vacuum polarization amplitude and gluon-coloron mixing amplitude. 
Particle notation is as defined in Fig. [S] 



F. Full Propagator Correction 



We have just seen that, due to the pinch technique, the coloron-coloron and gluon-coloron VPAs receive an addi- 
tional contribution from the pinched non-Abelian vertex corrections. Combining the VPAs, the UV divergences can 
be removed by two wavefunction renormalization counterterms (which arise from renormalizing the gauge eigenstates 
Gi/j^ and 6^2/1 ) and one mass counterterm (which arises from renormalizing the vacuum expectation value /), in the 
usual way. In the MS scheme we obtain 



\4ttJ 



(G) f 
Jo 



Ucc{q^)+Ucc{q'') 
+4cot^(2i9c)C2(G)^ dx' 



dx< 



^GC 



^GC 

1 j (^3 + 4x(l - x)^E + (1 - 2x 



ij 2(^3 + 4x(l + 2(1 -2x)^ 

-1 ) 2(2-x)E + 3-5x 



^cc 



i)Nf f 
Jo 



lJ^ + 1 

dx 



x(l — 9 
^ 



^cc 



-E 



Mi 



A, 



2x{l-x) 



E (f 



and 



.47r 



nGc(g') + nGc(^')^ 

= 2cot(2i9c)C2(G)^ dx\^ 



(57) 



(58) 



((£)'-) ki--'"-^-') 



Ml 



E^{\-2xf(^ 



dx 



A. 



1 



2x(\-x) E . 



The overall UV-finite propagator correction to the qq ^ C amplitude can be found by insering these expressions in 

(59) 



ir ^^pinched — l^J^^qq^c 9 "1" ^''^^qq^C 3 



Letting rrig ^ 0, we find that P + Ppinched becomes IR divergent, with the divergence arising from H-'cc Setting 



rria 



and e < gives 



^P + ^Ppinched = /C2(G) 

47r 



iM 



(0) 

qq^C 



finite . 



(60) 



We have seen that the pinched diagrams contribute to the full propagators of the gluon-coloron system. This 
might seem in confiict with the expectation that the mass poles should be a property of freely propagating particles, 
and should not depend on any initial and/or final state. However, when we sum the Dyson series to obtain the 
full propagator, the pinched diagrams always appear as an overall pre-factor, as pictorially shown in Fig. 11 This 
has an overall effect on the full propagators, which depend on the initial and final states, but has no effect on the 
propagator poles. Thus when we compute physical masses, we can do so by employing the "true" propagators in the 
computation, without the contribution from the pinched diagrams. 
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FIG. 11. The relevant contributions to the coloron Dyson series; as before, the zigzag hnes represent colorons. The first row 
is the sum of the coloron VPA diagrams in the propagator, while the second row represents the sum of the VPA diagrams 
on top of the one-loop contribution from the "pinched" vertex correction (the double curly line illustrates generically all the 
allowed gauge bosons in the original non-Abelian vertices). The overall pinched amplitude factors out, and has no effect on 
the coloron pole mass. 



G. Cross Section at One Loop 

Adding up the tree- level contribution and the NLO contributions from iQ + iVAbeiian + ^^non-Abeiian, and iP + 
^^pinched , givcs 3, qq ^ C amplitude of the form 



iM 



qq^C 



iMf^^C + + ^VAbelian + ^Kon-Abelian + + ^^^pinched 



^-^qq^C + 4^ 



(r iM 



qq^C ^ ^ ^-^^^qq^C ^ ^ ^-^^^qq^C I ' 



(61) 



where expressions for the real parts of T, T', and T" are given below. Averaging the squared amplitude over initial 
spins and colors, summing over final polarization states, and integrating over the phase space, gives the NLO result 
of the form 



.(0) 



f virt — '^qq'-^C + qq^C 



a,A{rl + rl) 



5{1-X) 



l.g Rer.^ 



4 I 4 



^ Re T' + "^^^^ Re T" 



' R 



■ (62) 



At 5 = it is possible to integrate over the Feynman parameter space in the expressions for iQ + ^VAbeiian + 
^Kon-Abeiian, and iP + iPpinched- As wc havc sccu, the UV infinities cancel in iQ + iVAbeiian + ^Kon-Abeiian and 
are absorbed by propagator conterterms in iP + iPpinched- Thus for rrig ^ the overall amplitude is finite. Taking 
the rrig limit leads to IR divergences in Re T, which are parametrized by taking e < 0. For small and negative 
values of e we obtain 



ReT=ri?iVr(l + e) 



+ 



Vm2 

77 Ttt 



2 3 „ 47r2 
8+ — 



C2{r) + 



^ 61 Stt 

- E-i ^ 

9 2V3 



29 , Ml 

log — — 

48 16\/3 16 ^ 



Re T' 
Re T" 



11 27r2 



C2{r) + 



4cot2(26'c)C2(G) 
57r 



27r2 

2\/3 3 



C2{G) , 



95 
~9 



iV^TT 43 , Ma 
^-12^°^^ 



2cot(26lc)C2(G) 



7r2 8 , M2 

^ - o log ^ 

3 3 fj,'' 



C2{G) 



1 I, Ml 
- + - log ^ 
9 6 /Lt^ 



irl + r%)Nf , 



1 , 

3l°^# 



(r-L + rR)Nf 



(63) 
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{a) (b) (c) 

FIG. 12. Diagrams contributing to the real emission processes. A gluon field is, as usual, represented by a coiling line; a 
coloron field is represented by a zigzag line, (a) Gluon emission, (b) Quark emission, (c) Antiquark emission. 



In the next section we will compute the corrections to the tree-level cross section due to the emission of soft and 
collinear gluons. We will show that the real emission cross section has IR divergences which exactly cancel the IR 
divergences contained in (Jvirt (Eq. (|62|), leading to a total cross-section free of both UV and IR divergences. 



V. NLO COLORON PRODUCTION: REAL CORRECTIONS 



The real emission corrections, at NLO, are given by the diagrams of Fig 12 We first consider the diagrams with 
real emission of a gluon, shown in Fig [12] (a). The squared amplitude, averaged over initial colors and spins, and 
summed over final colors and polarizations, is found to be, in d = 2{2 — e) dimensions. 



\M 



qq^gCl 



C2(r)fff(ri + 4) ^,^ 
dim(r) 



uj{l — uj) 



C2{r) + C2{G) 



where 



1 — cos i 



(64) 



(65) 



for the real gluon emission is 

where the integral is over the 2-body Lorentz-invariant phase space in parton CM. In d = 4 — 2e dimensions. 



1 1-x 



Sir r(l - e) 
This leads to the partonic cross section 



21-^ pi 



duj [uj (1 — uj)]' 



^(1) 



o^s{rl^rl)A /4Vy r(l - e) 



271 J r(l - 2e) 



2xHl±2(^ ^ 2x^(i + x + x') 



3 (1-x) 



l+2e 



(67) 



(68) 
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Now X is no longer constrained to be equal to one. Instead we must have x ^ oi" else no on-shell coloron can be 
produced. The term proportional to C2(r) features a collinear singularity, parametrized by e, and a soft singularity, 
parametrized by 1 — x- The term proportional to €2(0) only features a soft singularity. The integral over x in 
Eq. (68) is finite for e < 0, in spite of the singularity of the integrands. For small and negative values of e we can 
rewrite the x-dependence as follows 



x^(i + x') 

(l_^)l+2e 



1 1 + 

e (1-X)+ 

3 An A ^ 1 + X + x' 



2(1 + X') 



log(l - x) 
1-X 



1-X 



logX 



(1 - x)^+2^ 

where, as conventional, the distributions are defined by 



dx[fix)-fil)] 



log(l - X) 
1-X 



(69) 



(70) 



The coefficients of the delta functions are found by integrating both sides of the equations. The partonic cross section 
becomes 



(1) 

qq^gC 



asA{r 



27r 



5{i-x) R^R' 



(71) 



where, using Eq. (|50|), and expanding for small values of e, 



R 



r(l-e) 



C2(r) 



R' = -2 



£ - log ^ 



r(l-2e) 

P,^,(x) + C2(r) 



C2(G) 



1 



4(1 +X') 



2, / log(l - x) 



1-X 



1-X 



logX 



C2(G) 



2 1 + X + 

3 (1-X)+ 



(72) 



In the second equation Pq^q{x) is the Altarelli-Parisi splitting function for an on-shell quark to evolve into a virtual 
quark and a real gluon: 



q^q 



(X) = C2{r) 



(73) 



Adding together ^virt, given by Eqs. (62) and (63), and cr]^^^g(j^ given by Eq. (jrij), shows that the IR divergences 
proportional to 5(1 — x) cancel. There is still a collmear singularity in R' ^ proportional to the Altarelli-Parisi evolution 
Pq^q{x)- This singularity arises from integrating over all collinear initial-state gluons. As we will see in the next 
section, these collinear IR divergences will be absorbed through renormalization of the PDFs. 

The real quark and antiquark emission diagrams are shown in Fig. 12 (b) and (c), respectively. The corresponding 
summed-averaged squared amplitudes are 



X C2{r) + C2{G) 



C2{r)gt{rl+r%) 2. 



dim(G) 

(i-x)(i-^) 

(l-(l-X)(l-c^))^ 



(74) 



2(e + x) + 



1 



2x(l-x) 



(l-X)(l-c^) 



+ (l-e)(l-x)(l-w) 



and 



C2{r) + C2{G) 



C2{r)gi{rl+rl) 
dim(G) 

(1-x)'^ 



(75) 



(I-(I-X)c.)^ 



2(e + x) + 



1 



2x(l-x) 



(1-x) 



+ {1 - e){l - x) uj 



where dim(G) = 8 is the dimension of the adjoint representation. Note that the amplitudes for quark and antiquark 
emission are related by crossing, i.e. uj ^ (1 — uj). The integration over the two-body Lorentz-invariant phase space 
proceeds as in the gluon emission case, yielding 



(1) 

qg^qC 



cr. 



(1) ^ o^sAjrl + r\) Os 

qg^qC ~ s 2tT 



R" 



(76) 
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where 



R" 



dim(r) J 
dim(G) I 



C2(r) + C2{G) 



(l-x)(2 + x + 2x') 



X 



2(1 + x) logx 



(77) 



Here Pg^q{x) is the Altarelli-Parisi splitting function for an on-shell gluon to evolve to a virtual-real quark pair, 



dim(G) 



C.(r).dim(r) ^ _ 

where C2{r) • dim(r)/dim(G) = 1/2. There is no soft singularity in a^^^ 



'qg^qC 



a 



(1) 

qg^qC 



(78) 



, only a collinear singularity 



proportional to the Altarelli-Parisi evolution Pg^q{x)- As noted above regarding a^^i^^^^, this singularity will be 
canceled by renormalization of the PDFs when we compute the total hadronic cross section. 



VI. NLO CROSS SECTION 



Our calculations in the previous sections have produced all of the relevant partonic cross sections at NLO and 
demonstrated them to be both UV and IR finite. Note that the gg ^ C process vanishes at tree level ^ and 
the one-loop contributions are small, less than of order 0.1% of the g'^-initiated leading order contribution [35]; we 
therefore do not include it in this work. 

The full NLO cross section for coloron production at the LHC is 



NLO 



J dxi J (ix2 i 



C^l)fl{x2) + fl{xi)f°{x2) 



JO) 



q^gcj 



E 



q L 



a 



(1) 

qg^qC 



(79) 



where the partonic cross-sections a are given in Eqs. (20), (62), (71), and (76), and where the superscript "0" in 
the PDFs will be clear in a moment. We saw that all IR divergences contained in a cancel, except for a couple of 
collinear singularities proportional to Altarelli-Parisi evolutions. Such singularities arise because we integrated over 
all collinear quarks and gluons, even those which we should have included in the PDFs. Therefore, the corresponding 
IR singularities are absorbed by renormalizing the bare PDFs in Eq. (79). In the MS scheme. 



Mx, A^.) = /f (-) - 1^ ( ' - 7 + log(4.) - log ^) / f E (^) 



(X) 



(80) 



where z, j = g,^, and /uf is the factorization scale. Exchanging the bare PDFs for the renormalized ones replaces E 
with log in Eqs. (72) and (77). The hadronic cross section becomes 



NLO _(^s A Hi{Oc) 



f dxi f dx2 \ sr^ 
J Xi J X2 1 4^ 



fq{Xi,llF)fq{x2,llF) + f q{Xi, IIf) f q{x2, I^f) 
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E 



fq{Xi,llF)fg{X2,HF) + fg{Xi,llF)fq{.X2,tXF) + {fq ^ fq) 



27r 



(81) 



where the function Hi{Oc) is defined below, in Eq. (85), A is defined in Eq. (21), and the partonic CM energy s has 
been traded for the hadronic one, as in Eq. (25). Notice that since the integrand is now finite, we can ignore the 
1 — e factor in A. The functions T^^{x) and T^^{x) are 



c 



T^%X)=2\0g^.q^q 



^F 



(x) + ^,(x) , 



^«^(x) = log^P,^,(x)+D,(x) , 



(82) 
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where 
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(83) 



and 
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log^ COt(2^,)ff4(^c) 



(84) 



The functions Hi{6c) are determined by the chiral couphngs of the quarks to the colorons (which depend on the 
charges of the quarks under the fuh SU{3)ic x SU{3)2c symmetry): 



( 2 tan^ Oc 



vl =rR = - tan ( 



HsiOc) 



tan^ (9c + cot^ i9c vl ^ vr , i^2(<9c) = <( 

, 2 cot^ 6c = rR = cot Oc 
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tan"^ Oc + cot"^ Oc 



(2 rL = vr = - tdiii Oc 

2(l + cos(4i9c)) 
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L cot^ 
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3 + cos(4(9c)) 

— cot Oc 

sin(4^,) 
3 + cos(4i9c) 

, tan 6>c 



tl + tr 

rL = rR = cot Oc 
Vl =rR = - tan Oc 



VL^VR . (85) 



Vl = rR = cot ^ 



At NLO the fi dependence is removed by trading the MS couphngs gis and g2s^ or gs and Oc^ for the corresponding 
running couphngs. Since Oc is a free parameter, we simply set /i = Mc, and express the cross section as a function of 
the MS couphngs. At the same time, the NLO jj^f dependence weakens once the renormahzed PDFs are employed, 
as a in Eq. (79) is independent of /ii? to this order in perturbation theory. 



From these results we may also compute the transverse momentum distribution of the produced coloron, which is 
given by 
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FIG. 13. Behavior of the Q function defined in Eq. ( |84| , for fi — Mc'- this gives the contribution from the virtual corrections 
to the NLO cross section for coloron production. The upper curve is for the tl ^ tr scenario, whereas the almost identical 
lower curves are for vl = vr^ — tan^c and tl = tr^ cot^c- Note that Q, and therefore the NLO corrections, become very 
ge when sin^ Oc is either too small or too large. 
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FIG. 14. Dependence of LO and NLO cross sections at the LHC (y/s = 7 TeV), as a function of factorization scale /jf for 
Mc — 2.0 TeV, sin^ ^cU=2.0TeV = 0.25, and the three possible flavor- universal scenarios for the quark charge assignments. As 
expected, the NLO cross section has a much weaker (formally, two-loop) residual scale-dependence. 



and uj (Eq. (65)) is given by 



1- Wl 



(89) 

Note that this is the leading ort/er prediction for da/dpr^ and therefore this distribution is strongly /ii?- dependent. 



VII. DISCUSSION 



We now illustrat^ our results for the NLO coloron production cross-section in Figs. 13 - 16 In each figure we 
consider the three possible flavor-universal scenarios for quark charge assignment: vl = vr = — tan 6>c, vl / ^i?, and 
vl = rR = coiOc All of the plots refer to coloron production at the LHC with y/s = 7 TeV. 

Notice that the perturbative expansion is only meaningful as long as sin^c is neither too close to zero (where 
92s ^ gis) nor too close to one (where gis ^ ^2s)- This is clear from Fig. 13, in which we plot the quantity Q defined 
in Eq. (84), for /i = Mc: the contribution from the virtual corrections to the NLO cross section. The upper curve is 
for the tl scenario, whereas the almost identical lower curves are for = vr = — tan6>c7 and vl = tr = cot 6>c- 
For sin^ Oc ^ 0.05 and sin^ Oc ^ 0.95 the virtual corrections become large, and the perturbative expansion in as 
breaks down. Since 0.118 at the Z pole, these boundaries correspond to g2s ^ 2.7 and gis ^ 2.7, respectively. 

we plot the fiF dependence of the LO and NLO production cross sections of a 2.0 TeV coloron (with 



In Fig. 



14 



TeV=0.25). The scale-dependence of the LO cross section is of order 30% while, as expected, the NLO 



cross section has a much weaker scale dependence, only of the order of 2% percent. 



^ For the purposes of illustration we use the Mathematica package for CTEQ5 36 to evaluate the relevant parton distribution functions. 
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FIG. 15. NLO cross section times branching-ratio to quarks for on-shell coloron production at the LHC (y^ = 7 TeV), 
corrected for acceptance as described in the text. We consider the three possible flavor-universal scenarios for the quark 
charge assignments, take the renormalization scale /i to be equal to Mc, and plot a for sw? Oc\iJi=Mc — 0-05 (dashed), 0.25 
(dot-dashed), and 0.5 (dotted). We plot these cross sections for fip ranging from Mc/'^ to 2Mc and, reflecting the weak 
dependence of the NLO cross section on the factorization scale, the resulting bands for each sin^ Oc are very narrow. To give 
a sense of current experimental reach, we plot the CMS [22] upper limit (solid line) on the cross-section times dijet branching 
ratio for a narrow resonance. Note that the axigluon [T] corresponds to the middle tl ^ tr plot with sin^ ^ = 0.5 - and hence 
a narrow axigluon resonance is constrained to have a mass of order 2.6 TeV or higher. 
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FIG. 16. "iT-factor", the ratio of the NLO to LO cross section for coloron production at the LHC {^/s = 7 TeV), plotted as 
a function of Mc for sin^ = 0.05 (dashed), 0.25 (dot-dashed) and 0.50 (dotted), /jf = Mc, and the three different quark 
charge assignments. 

In Fig. [15] we plot the cross section times branching ratio to quark jets as a function of Mc, allowing /j^f to vary 
from Mc/2 to 2Mc. Here, in order to compare to the experimental results of [22 (shown as the solid line in the 
figures), we correct for the acceptance of the detector by multiplying our partonic- level NLO production cross section 
by the factor 



R = 



{a{pp C) 



axigluon 



axigluon 



(90) 



In this expression, {(j{pp C) 'B'A)a:^uon CMS (LO) prediction for axigluon production cross section, times 

dijet branching ratio, times acceptance^ reported in [22 , and a^^{pp C) axigluon is the leading order cross section 
in Eq. (23) in the case of an axigluon (J7e. tl = —tr = 1), assuming the branching ratio to quarks B{C qq) = I p] 
The three sets of thin bands correspond to sin^ ^c|/i=Mc = O-O^ (dashed), 0.25 (dot-dashed), and 0.5 (dotted). Here, 
the weak residual /j^f dependence is shown by the narrowness of the bands. To give a sense of current experimental 
reach, we also show the 1 fb~^ CMS upper bounds on the cross-section times di-jet branching ratio for a narrow 
resonance JT^ 



Note that the bound on the axigluon [1 corresponds to the vl ^ tr plot with sin = 0.5 - and 
hence a narrow axigluon resonance is constrained to have a mass of order 2.6 TeV or higher. The enhancement of 
the axigluon cross section at NLO is responsible for the increase in the bound from of order 2.5 TeV as reported in 

Next, we compute the "i^-factor" for coloron production. 

K{Mc,smOc\^=Mc.f^F = Mc)- ^ 



fc,Sm (ic\fi=Mcil^F 



Mc) 



a^o^Mc,smOc\^=Mc.l^F = Mc) 



(91) 



The CMS acceptance for isotropic decays is of order 0.6, independent of resonance mass [22| . 

It is worth noting that there are examples of models with colorons which do not decay primarily to dijets, e.g. 
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shown in Fig. 16 for sin^ = 0-05 (dashed), 0.25 (dot-dashed) and 0.50 (dotted). Again, we see that the NLO 
corrections are of order 30%. In Appendix |b| we report the numerical values of the iC-factors corresponding to Fig. 
16, as well as those corresponding to the ATLAS KK-gluon search reported in [38] , 




FIG. 17. Fraction of colorons produced with a pr greater than pTmin, as a function of pTmin- The curves are for Mc — 1.2 
(highest, thin blue line), 2.0 (middle, medium purple line), and 3.0 TeV (lowest, thick green line), for the vectorial case 
tl — tr — — tan^c and sin^ 6c — 0.05. Note that of order 30% of the colorons in this mass range are produced with pT > 200 
GeV. As denoted by the red shaded region, below a pr of 200 GeV the corrections become larger than 30%, terms proportional 
to \og{Mc / p'^^^n) become large, and this fixed-order calculation becomes unreliable. 



At leading order, the coloron is produced with zero transverse momentum. We may use our results to compute 



the pt spectrum in coloron production to leading non-trivial order from Eq. (86). Using these formulae, we may 
compute the fraction of colorons produced above a momentum pTmin 

cr^^<^(Mc, sm i9c I ;.=Mc,Mf = Mc) dpT 

where prmax is the kinematic maximum transverse momentum (which depends on the coloron mass). For illustration, 
we plot this fraction for vectorial colorons (r^ = tr = — tan 6>c, with sin^ Oc = 0.05) with masses of 1.2, 2.0, and 3.0 



TeV in Fig. 17 Note that of order 30% of the colorons in this model and mass range are produced with pt > 200 
GeV. Below a pr of 200 GeV the corrections become larger than 30%, terms proportional to \og{M^ / p'^^-^) become 
large, and this fixed-order calculation becomes unreliable. 

In conclusion, we have reported the first complete calculation of QCD corrections to the production of a massive 
color-octet vector boson. Our next-to-leading-order calculation includes both virtual corrections as well as corrections 
arising from the emission of gluons and light quarks, and we have demonstrated the reduction in factorization-scale 
dependence relative to the leading-order approximation used in previous hadron collider studies. In particular, 
we have shown that the QCD NLO corrections to coloron production are as large as 30%, and that the residual 
factorization scale-dependence is reduced to of order 2%. We have also calculated the i^-factor and the pt spectrum 
for coloron production, since these are valuable for comparison with experiment. Our computation applies directly to 
the production of the massive color-octet vector bosons in axigluon, topcolor, and coloron models, and approximately 
to the production of KK gluons in extra-dimensional models or colored technivector mesons in technicolor models. 
We look forward to future results from the LHC, and the possible discovery of colorons. 
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Appendix A: Feynman Rulej^ 



The Feynman rules for the trihnear and quartic vertices are shown in Figs. 18 and [T9j respectively The coloron 
is represented by a zigzag line, the coloron ghost by a sequence of small circles, and the eaten Goldstone bosons by 
dashed lines. All other particles are denoted as in QCD standard notation. 
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FIG. 18. Feynman rules for the trilinear vertices. In each diagram the momenta are toward the vertex. A gluon field is, 
as usual, represented by a coiling line; a coloron field is represented by a zigzag line. The coloron ghost is represented by a 
sequence of filled circles, and the eaten Goldstone bosons are represented by dashed lines. 



Appendix B: Numerical Values of the i^- Fact or 



The numerical values of the i^-fa c tors fo r var ious va lues of the coloron mass and the three patterns of coloron 
coupling are shown in Figs. 2Q(a)[ |2Q(b)[ and 2Q(c)[ Finally, the values of the K-factor corresponding to the 
KK-gluons of [12], corresponding to the experimental search reported in are shown in Fig. 21 
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^ The Feynman rules discussed here are equivalent to those in ^24^ , aside from those for the triple-coloron vertex which is not specified 
in that reference. 
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FIG. 19. Feynman rules for the quartic vertices. A gluon field is, as usual, represented by a coiling line; a coloron field is 
represented by a zigzag line. 
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FIG. 20. K-factors for colorons of various masses and couplings. The classic "axigluon" 1 corresponds to ^ vr and 
sin^ Oc = 0.50. 
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FIG. 21. X- factors for KK-gluons of various masses considered in [38^. This calculation is based on the theoretical framework 
of [12], with the KK-gluon coupling (specified in the column heading) varying between — 0.20^s and —0.40^^. 
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